Jahn- Teller stabilization of a "polar" metal oxide surface: Fe 3 O4(001) 
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Using ab initio thermodynamics we compile a phase diagram for the surface of Fe3O4(001) as a 
function of temperature and oxygen pressures. A hitherto ignored polar termination with octahedral 
iron and oxygen forming a wave-like structure along the [110] -direct ion is identified as the lowest 
energy configuration over a broad range of oxygen gas-phase conditions. This novel geometry is con- 
firmed in a x-ray diffraction analysis. The stabilization of the Fe3O4(001)-surface goes together with 
dramatic changes in the electronic and magnetic properties, e.g., a halfmetal-to-metal transition. 

PACS numbers: 68.35.Md, 68.35.Bs, 68.47.Gh, 73.20.At, 75.70.Rf, 61.10.Nz 



The surface composition and structure of a multicom- 
ponent material is a fundamental property that has sig- 
nificant consequences for its surface reactivity, mechani- 
cal, and magnetic properties. Magnetite is not only im- 
portant in geophysics and mineralogy, but is also attract- 
ing increasing attention as a potential material for spin- 
tronic devices Q, Q], due to its half-metallic behavior, 
coupled with a high Curie temperature of 858 K. 

Magnetite crystallizes in the inverse spinel structure, 
where the oxygen atoms form a slightly distorted face 
centered cubic lattice and the iron atoms occupy tetra- 
hedral and octahedral interstitial sites. The stacking se- 
quence in the [001]-direction consists of A-layers contain- 
ing tetrahedral iron (FeA) and B-layers with oxygen and 
octahedral iron (Fee) atoms (cf. Fig. 

Two alternative approaches are widely used to assess 
possible terminations for metal oxide surfaces, an ionic 
scheme considering the electrostatic energy |2j and the 
autocompensation rule applied to covalent bonds 0] . Af- 
ter the scheme of Tasker Q, bulk truncations of the 
Fe3O4(001)-surface should be unstable because of a di- 
verging surface energy due to uncompensated dipole mo- 
ments. Such structures were, therefore, so far discarded. 
Both the ionic Q and the covalent Q approach rely on 
the information from the bulk only and suggest that sta- 
bilization of a polar surface can only be achieved through 
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FIG. 1: The inverse spinel structure of magnetite together 
with a top view of the two bulk truncations of Fe3O4(001) 
with an A- and B-layer, respectively. Oxygen atoms, FeB and 
FeA are marked by white, grey, and black circles. 



massive changes of the surface stoichiometry. 

For Fe 3 O 4 (001) a (a/2 X \/2)J?45 -reconstruction has 
been observed experimentally 0, IH El El Ell • How- 
ever, the surface structure and stoichiometry are still 
under debate. Applying the autocompensation rule Q 
to Fe 3 O 4 (001), two models for the (a/2 x \/2)i?45°- 
reconstruction have been proposed so far: an A-layer 
where half of the tetrahedral iron is missing 0,0>0and a 
B-layer with oxygen vacancies or hydroxyl groups 0, ^| . 
Several variations of the former model have been sug- 
gested: The x-ray photoelectron diffraction (XPD) anal- 
ysis of Chambers et al. and the low energy ion scatter- 
ing (LEIS) results of Mijiritskii et al. E| were interpreted 
as a half-filled A-layer with strong inward relaxations of 
the surface layers. Concerning the magnitude of the in- 
terlayer distances, there is a considerable discrepancy. In 
a molecular dynamics calculation with classical poten- 
tials - the only theoretical work so far - Rustad et al. fll| 
suggested a 0.5ML A-termination where the surface and 
half of the subsurface FeA relax in octahedral positions in 
the surface B-layer. Also non-autocompensated termina- 
tions with a charge ordered B-layer were reported from 
scanning tunneling microscopy (STM)-measurements us- 
ing a ferro- |l2| or an antiferromagnetic tip |l3| . However, 
the observed magnetic contrast alone cannot explain the 
LEED- and XRD-superstructure intensities and does not 
provide a clue on the stabilization mechanism of the sur- 
face. 

In this Letter we present the results of a systematic in- 
vestigation of the composition, structure and properties 
of the Fe3O4(001)-surface based on density-functional 
theory (DFT) calculations. Using ab initio atomistic 
thermodynamics 0, 0] , we predict that a "polar" ter- 
mination is the lowest energy configuration over the en- 
tire range of accessible oxygen pressures. The stabi- 
lization of the surface involves a fundamentally differ- 
ent mechanism, which has not been considered so far: 
While most of the previous studies proposed an ordering 
of surface vacancies as the origin of the surface recon- 
struction, here it is explained as a Jahn- Teller distortion 
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of the surface atoms forming a wave-like pattern along 
the [110]-direction. A x-ray diffraction (XRD)-analysis 
supports this theoretically predicted model. 

We use the full-potential augmented plane waves (FP- 
APW) method in the WIEN2k-implementation [lffl and 
the generalized gradient approximation (GGA) [ijj . The 
Fe3O4(001)-surface is modeled by a symmetric slab con- 
taining five B- and four to six A-layers depending on 
the structural model. The vacuum between the repeated 
slabs amounts to 10 A. We have ensured convergence 
with respect to the thickness of the slab by repeating 
the calculations for the most stable configurations with 
a thicker slab with an additional A- and B-layer on both 
sides (approximately 100 atoms in the unit cell as op- 
posed to 70 in the original unit cell). The surface free 
energies for both supercells are within ±3meV/A 2 equal. 
The lateral parameter of the supercell is set to the GGA 
bulk lattice constant, 8.42 A, which is in good agreement 
with the experimental value of 8.394 A. 

The convergence parameters for the mixed APW+lo 
and LAPW basis set are: = 1.90 bohr, R^ T = 

1.40 bohr, inside the muffin tins (MTs) wave functions 
expansion in spherical harmonics up to l^ax = 10 ano - 
non-spherical contributions to the electron density and 
potential up to = 4. The energy cutoff for the plane 
wave representation in the interstitial is = 19 Ry for 
the wave functions and — 196 Ry for the potential. 

With these cutoff parameters a convergence of energy 
differences better than 1 mRy is achieved. Results for the 
(\/2 x v / 2)i?45°-unit cell are obtained with 4 fc||-points in 
the irreducible part of the Brillouin zone. 

The configurations we have considered include (i) a 
0.5 ML of tetrahedral iron above the B-layer H 0, H, 
(ii) or relaxed in the B-layer; (Hi) a modified 0.5 ML A- 
termination where the surface and half of the subsurface 
FeA occupies octahedral sites in the B-layer ^l|. Addi- 
tionally, we have considered two geometries with one oxy- 
gen vacancy per unit cell: either as a nearest neighbor to 
a tetrahedral iron or above an octahedral iron. Together 
with the non-autocompensated A- and B-terminations 
shown in FigE]they were used as a starting point for the 
DFT-calculations, where we performed a full structural 
optimization of the atomic positions in the outer two AB- 
double-layers with damped Newton dynamics 0] . 

The calculated surface phase diagram for all studied 
terminations is compiled in Fig. [5] The oxygen chem- 
ical potential is converted into oxygen pressures for a 
fixed temperature of T = 900 K, which corresponds 
to a typical annealing temperature in experiment. The 
range of allowed oxygen chemical potentials lies between 
the "oxygen-poor" -limit, which corresponds to the case 
where the oxide would decompose in iron crystallites and 
oxygen gas and the "oxygen-rich" -limit, defined by the 
chemical potential of oxygen in the oxygen molecule, 
set as zero reference. In Fig. [5] the surface energy is 
given with respect to bulk magnetite. We note that at 



B-layer modified 




H [eV] 

FIG. 2: Calculated surface free energy 7(T, p) as a function of 
the chemical potential of oxygen (bottom x-axis) for all stud- 
ied terminations. In the top x-axis the fio is converted into 
a pressure scale at T = 900 K. The vertical lines mark the 
oxygen-poor and oxygen-rich limits of the oxygen chemical 
potential, with 1 /2_Z?o 2 taI used as zero reference. The dashed 
region marks the range of pressures that were used during 
sample preparation in the experiment. The B-termination 
with bulk atomic positions and modified positions with a 
(a/2 x v / 2)-R45°-periodicity are marked as ideal and modi- 
fied, respectively. The B-layer with oxygen vacancies above 
an octahedral iron and next to a tetrahedral iron are denoted 
by (1) and (2), respectively. Top and side view of the mod- 
ified B-layer is given at the top. Surface oxygens with and 
without a FeA neighbor are denoted by 02(S) and Ol(S). For 
the color code see Fig. 

T = 900 K and po 2 > 10~ 6 mbar hematite becomes the 
thermodynamically stable phase 0|. Taking into ac- 
count the calculated enthalpy of formation for magnetite 
(11.32 eV) and hematite (8.024 eV) j^, the stability 
range of Fe304 extends up to fio = — 1-5 eV. 

Autocompensated terminations, e.g., the ones contain- 
ing 0.5 ML FeA m the surface layer, are independent of 
the oxygen pressure, while the effect of the oxygen en- 
vironment is most distinct for the two bulk truncations 
with an A- or a B-layer. In line with electrostatic con- 
siderations, the A-tcrmination is unstable and this trend 
is substantially enhanced for oxygen-rich conditions. A 
most striking result is that the other bulk truncated sur- 
face, the B-termination, which has been so far discarded 
as a non-compensated surface, turns out to be the most 
stable configuration over a broad range of oxygen pres- 
sures even in the ideal (unrelaxed) geometry. For oxygen- 
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FIG. 3: Results of the XRD-analysis. The experimental data 
is marked by squares. Solid and dashed lines correspond to 
fits assuming the modified B-termination and 0.5 ML FeA- 
termination, respectively. 



poor conditions it competes with a 0.5 ML tetrahedral Fe 
relaxed in the B-layer and with a B-termination with O- 
vacancies lying above an octahedral iron in the subsurface 
B-layer. An asymmetric force on the surface oxygens in 
the unreconstructed ideal B-termination entails a wave- 
like-distortion along the Fee-row. Hence, the origin of the 
symmetry reduction and the observed (y/2 x y/2)RA5°- 
pattern is a Jahn- Teller-distortion [2ll |22| . Compared 
to the ideal B-termination, the formation of this modi- 
fied B-termination is connected with an energy gain of 
approximately 20 meV/A 2 and a reduction of the work 
function of 0.46 eV from 5.78 eV to 5.32 eV. We note that 
a reduction of <t> indicates a decrease of the dipole mo- 
ment induced on the surface, i.e., a decrease in polarity. 
The calculated work function is in good agreement with 
the value measured in photoemission experiments 23]. 

A top and side view of the novel surface geometry, is 
given at the top of Fig. [21 The wave-like-pattern along 
the [110]-direction emerges from a strong lateral relax- 
ation of 0.30 A of one of the two oxygen atoms with a 
missing FeA-neighbor (marked as 01 (S) in Fig. EJ to- 
wards the Fee-row and a subsequent alternating shift 
of Fee-pairs by ±0.09 A perpendicular to the Fee-row. 
The oxygen atoms with a FeA-neighbor (02(S)) relax by 
0.14 A towards and by 0.03 A parallel to the Fee-row. 
Due to the lateral relaxation of oxygen the Fee-oxygen 
distance is reduced from 2.06 A in the bulk to 1.96 A on 
the surface and are thus closer to the bulk values for the 
tetrahedral (^o-Foa = 1-89 A) coordination. Perpen- 
dicular to the surface, all surface atoms relax inwards 
except for 01 (S), which relaxes outwards. The atoms in 
the subsurface A- and B-layer relax outwards. A detailed 
structural analysis will be presented elsewhere 

X-ray measurements were performed on a natural crys- 



tal [25| at the beam line DW12 at LURE/Orsay with a 
wave length of 0.826 A. In total, 220 symmetrically inde- 
pendent reflections were measured at room temperature 
in 5 crystal truncation rods (CTR) and 6 superstructure 
rods. The cleanliness of the sample and the reproducibil- 
ity of the (a/2 x \/2)i?45 -LEED pattern were verified 
by a repeated sputtering and annealing up to 850 K in 
UHV and oxygen pressures up to 5 x 10~ 6 mbar. Fur- 
thermore, the LEED I/V curves of the natural and an 
artificial crystal are nearly identical |2fiT |. 

To minimize the effect of surface inhomogeneities only 
the superstructure reflections were considered in the anal- 
ysis. The experimental data and the calculated curves 
for the two most promising terminations, i.e., the half- 
filled A-layer of tetrahedral iron and the modified B- 
termination are displayed in Fig. [3J The corresponding 
R-factors obtained after a structural refinement with the 
DFT-geometries are 0.24 and 0.19, respectively, which in- 
dicates a preference for the modified B-termination. Both 
models vary only by the additional tetrahedral iron in 
the half-filled A-layer. The relaxations of the remaining 
atoms are rather similar leading to an overall similarity in 
the intensity distribution along q z . Still the fit with the 
modified B-layer is noticeably better. The B-termination 
with oxygen vacancies, proposed by Stanka et al. || is 
incompatible with the XRD-data (R-factor worse than 
0.6). The wave-like structure along the [110]-direction 
observed in the STM-measurements of Stanka et al. 
is a further evidence for the modified B-tcrmination 

The stabilization of the Fe3O4(001)-surface and the 
unusual lattice distortions are connected with strong 
changes in the electronic and magnetic properties. In 
the bulk, magnetite shows a half-metallic behavior with 
a band gap in the majority spin channel of approximately 
0.5 eV and a 100% spin-polarization due to the i2g-states 
of Fee at Ep in the minority spin channel |28j . 

On the surface (cf. Fig. surface states in the ma- 
jority spin gap reduce the spin-polarization at the Fermi 
level to approximately 60%. Spin-resolved photoemission 
measurements support this result (23. Since electrons 
near the Fermi level are important for transport proper- 
ties, we examined closer their nature: the states at Ep 
are confined to the surface layer and are due to the hy- 
bridization between d X2 _ V2 -states of Fee and p^p^-states 
of 01 (S). Additionally, there is a noticeable contribution 
from the tetrahedral iron in the subsurface layer. In order 
to get more insight into the driving mechanisms towards 
the reconstruction of the surface we have plotted in Fig.Q] 
also the difference between the DOS of the B-termination 
with modified and ideal positions of the atoms. For the 
modified B-termination a band gap opens in the majority 
spin channel between —0.4 and —0.2 eV where states of 
d X2 -y 2 (d Z2 ) character are shifted to higher (lower) en- 
ergies. The latter cross the Fermi level and eventually 
become partially unoccupied. Vice versa, in the minor- 
ity spin channel, bands become occupied. Altogether this 
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analysis can be misleading. 
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FIG. 4: Total density of states (DOS) of Fe 3 O 4 (001) termi- 
nated by a modified B-layer for the majority (minority) spin 
channel in the upper (lower) panel is shown in grey. The 
difference of the DOS of the modified B-termination and the 
bulk-terminated B-layer is given with a black line. A Gauss 
broadening of 2a = 0.05 eV is used. 



reduces the energy of the system. 

This combination of a Jahn- Teller-distortion and a 
spin-flip leads to interesting magnetic properties. In 
the bulk magnetite is a ferrimagnet with the magnetic 
moments of tetrahedral (—3.43 /ib) and octahedral iron 
(3.50 ^b) oriented antiparallel to each other and a total 
magnetic moment of 4.0 per formula unit. In the mod- 
ified B-termination the magnetic moments of the subsur- 
face Fsa (— 3.45^b) is close to the bulk value. However, 
the magnetic moment of the surface FeB (3.00 /is) is sub- 
stantially reduced due to the strong relaxations. 

Additionally a substantial magnetic moment of ap- 
proximately 0.25 /ib is induced in the undercoordinated 
surface oxygen with a missing bond to tetrahedral iron, 
while the magnetic moments of the rest of the sur- 
face oxygen atoms are close to the bulk value, 0.07/ib- 
A similar magnetization was observed for the oxygen- 
terminated a-Fe 2 O 3 (0001)-surface |2(ij . 

Based on ah initio thermodynamics we identify a B- 
termination with octahedral iron and oxygen in the sur- 
face layer forming a wave-like structure along the [1 In- 
direction as the thermodynamically stable configuration 
over a broad range of oxygen pressures. Instead of 
ordering of surface defects the symmetry reduction is 
achieved by a Jahn- Teller distortion and is connected 
with strong changes of the electronic and magnetic prop- 
erties, e.g. a metalization at the surface. Experimen- 
tal evidence for the DFT-predicted geometry is given 
in STM-measurements 0, |2j| and a XRD-analysis per- 
formed with the atomic positions obtained with DFT. 

Theoretically, "polar" surfaces have been predicted for 
a number of metal oxide surfaces, e.g. Fe2O3(0001) [20| . 
RuO 2 (110) 0, ZnO(0001) 0, and PdO HJ. The 
modified B-termination of FeaO^OOl) is a further ex- 
ample that excluding such surfaces from the structural 
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